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Universality Class of Ferromagnetic Transition in 
Three-Dimensional Double-Exchange System 
- 0(N) Monte Carlo Study - 



O 

o 

< 

(N 



I 

o 
o 



in 
in 

o 
m 
o 



o 
o 



X 



Yukitoshi Motome and Nobuo FURUKAWA^ 

RIKEN (The Institute of Physical and Chemical Research), Wako 351-0198 
^Department of Physics, Aoyama Gakuin University, Sagamihara 229-8558 

(Received ) 

Curie temperature and exponents are studied for the three-dimensional double-exchange 
model. Applying the 0{N) Monte Carlo algorithm, we perform systematic finite-size scaling 
analyses on the data up to 20'^ sites. The obtained values of the critical exponents are consistent 
with those of the Heisenberg universality class, and clearly distinct from the mean-field values. 
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Rediscovery of colossal magnetoresistance (CMR) phe- 
nomena has revived much interests in manganese ox- 
ides.^' As sample quality and experimental precision 
are improved, many experiments have been done to ex- 
amine critical properties of the ferromagnetic transition 
in these compounds. Critical properties, especially crit- 
ical exponents which assign the universality class of the 
transition, are highly important to understand the na- 
ture of the transition as well as basic physics of the CMR 
phenomena. Estimates of critical exponents, however, 
are scattered between different experimental results and 
widely range from the Heisenberg-like ones to the mean- 
field-like ones.'^^ There still remains experimental con- 
troversy on this issue. 

Theoretical estimates of critical exponents in the 
double-exchange (DE) model, which is considered to be 
a relevant model in these compounds,'*^ would help to 
reconcile this experimental controversy. However, it is 
difficult to clarify critical properties of this transition be- 
cause the system is a strongly-correlated electron system; 
itinerant electrons interact with localized spins through 
a strong Hund's-rule coupling whose magnitude is much 
larger than the bandwidth of electrons. Near the finite- 
temperature transition, spin fluctuations are critically 
enhanced due to the interplay between itinerant electrons 
and localized spins. A sophisticated tool is necessary to 
handle this many-body problem. 

Numerical calculation which includes all the correla- 
tion effects is one of the most promising tools for this 
problem. Monte Carlo (MC) method has been applied 
for studying the critical properties of this system. How- 
ever, most of the MC studies suffered from finite-size ef- 
fects since they were limited to small-size systems due to 
the rapid increase of the computational cost as the sys- 
tem size increases. Systematic analysis on the finite-size 
effect is indispensable to examine the critical properties 
of this transition in detail. 

The authors studied this problem by developing a 
MC algorithm, which is called the polynomial expansion 
Monte Carlo (PEMC) method. 5' '^^ In this method, the 
computational cost has reduced to 0{N^) from 0(A^^) 
in the conventional MC method^-* {N is the system size). 



and the finite-size scaling was applied to results up to 
8^ sites. The problem was also studied by Alonso et 
a/.^' They applied the hybrid MC algorithm in which 
the computational cost is proportional to 0{N^) (they 
claim it can be reduced to 0{N) empirically), and es- 
timated the critical exponent from the data up to 12^ 
sites. In both studies, the estimates of the exponents 
are consistent with those in the universality class of the 
Heisenberg model with short-range interactions. They 
are, however, not enough to estimate the values precisely 
and to show a clear deviation from the mean-field val- 
ues. Further study supported by an improvement in the 
algorithm is strongly desired for this purpose. 

In this work, we give precise estimates of critical ex- 
ponents of the ferromagnetic transition in the three- 
dimensional (3D) DE system. We apply an improved 
PEMC algorithm in which the computational cost scales 
to 0{N),^'> and perform the systematic finite-size scaling 
for the data from 6^ to 20^ sites. 

The DE model considered here is defined by the Hamil- 
tonian"*) 



{ij).o- 



h.C.)~ JH^Cr,.5,;, (1) 



where the first term describes the nearest-neighbor hop- 
ping of itinerant electrons on 3D cubic lattices and the 
second term is for the Hund's-rule coupling between the 
electrons and the localized spins 5'. For simplicity, we 
take the limits of Jjj — > cxs and S — \S\ — > cxd (classical 
spins) in the following calculations. 

We study the ferromagnetic transition in model (1) 
by using an improved PEMC method. In this method, 
we introduce effective truncations for the vector-matrix 
product and the trace operation, which reduce the com- 
putational cost for one MC update to 0{N). Readers are 
referred to Ref. 9 for the details of this algorithm. This 
technique enables us to study up to 20'^ sites within a rea- 
sonable time on parallel computers. We choose periodic 
boundary conditions and have typically run 4000 MC 
samplings for measurements after thermalization with 
1000 MC steps. The doping concentration is fixed at 
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quarter filling, namely, 0.5 electron per site on average. 
The energy unit is the half-bandwidth of noninteracting 
case, W = 6t = 1. 

Figure 1 shows the temperature dependence of the 
magnetization m = [S{k = 0)/7V]^/^. Here S{k) is 
the spin structure factor which is defined as S{k) = 



Sj) exp (ifc • rij)/N, where the bracket denotes 



the thermal average for the grand canonical ensemble. 
The magnetization data are the thermodynamic-limit 
values which are estimated by the system-size extrap- 
olation of S{Q)/N as shown in the inset. We fit the 
magnetization data by assuming the scaling relation 
TO oc (Tc — T)^ (the gray curve in the figure) to esti- 
mate the value of Tc and the critical exponent (3. The 
fit gives Tc = 0.0226(2) and f3 = 0.36(1). (Errors in 
the last digit are shown in parentheses.) The estimate of 
P is consistent with the Heisenberg one (3 = 0.365 and 
definitely distinct from the mean-field one /3 — 0.5. 

We also examine the finite-size scaling plot for S{Q) 
to determine the critical exponents. The scaling form is 
given by S{0)L^~-^ — f(eL^/'^) which is derived based on 
the hyperscaling hypothesis. Here e = [T — Tc)/Tc, L = 
7V^/^ and / is the scaling function. We plot S'(0)L''~^ as 
a function of eL^/'^ , and optimize the values of {Tc,7'i, v) 
so that all the MC data for different L and T converge 
to a single curve. Figure 2 shows the optimized fit. 
This fit gives the estimates as Tc = 0.0222(-H5, -3), v = 
0.60(-fl3, -10) and rj = 0.07(+34, -7). By using the 
relation /3 = z/(l -|-7y)/2 with the obtained distribution of 
f and 77, we obtain the estimate /3 — 0.32(-f 10, —6). Both 
estimates of Tc and /3 are consistent with the estimates 
in Fig. 1. We confirmed that the mean- field exponents 
f = 1/2 and = do not give a satisfactory convergence 
in the scaling fit. 

These results indicate that the ferromagnetic transi- 
tion in the DE model (1) belongs to the universality class 
of the Heisenberg spin system with short-range interac- 
tions rather than to that of the mean-field solution. This 
suggests that the critical exponents in the Heisenberg 
universality class should be observed in manganites when 
the DE mechanism plays a key role in the transition as 
widely believed. Other additional elements which are not 
included in the simple model (1), such as the electron- 
lattice coupling, may not affect the universality class be- 
cause the fluctuations through the additional degrees of 
freedom will make correlations more short-ranged. 

To summarize, we have investigated the critical ex- 
ponents in the ferromagnetic transition in the three- 
dimensional double-exchange model by using the im- 
proved Monte Carlo method. The finite-size scaling anal- 
ysis up to 20^ sites indicates that the transition be- 
longs to the universality class of the Heisenberg spin 
system with short-range interaction and is distinct from 
the mean- field one. In manganese oxides, we expect 
this Heisenberg universality class to be observed if the 



double-exchange mechanism plays a major role in the 
transition. Our results give a theoretical background to 
reconcile the experimental controversy on the estimates 
of the critical exponents. 
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Fig. 1. Temperature dependence of the magnetization on the ther- 
modynamic limit. The gray curve is the least-squares- fit to 
m oc (T — Tc)^ ■ Inset: System-size extrapolation of the spin 
structure factor. The data are for T = 0.013,- ■ -,0.0225 from 
top to bottom. 
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Fig. 2. The best-fit result of the finite-size scaling plot for the 
Monte Carlo data. 



